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SUMMARY 


A relatively simple and systematic method for forming the products of multiple 
Fourier series using tensor-like operations is demonstrated. This symbolic multi- 
plication can be performed for any arbitrary number of series, and only the coeffi- 
cients of one series need to be known. The application of this methodology to the 
transformation of a set of linear differential equations with periodic coefficients 
from a rotating coordinate system ;.o a nonrotating system is also demonstrated. It 
is shown that using Fourier operations to perform this transformation make it easily 
understood, simple to apply, and generally applicable. 


INTRODUCTION 


The harmonic balance method is a useful tool for obtaining the forced response 
of systems governed by differential equations with periodic coefficients. However, 
for a set of linear equations, the rno.^-ual application of this method tends to become 
quite unwieldy as the number of harmonics, or the number of degrees of freedom, 
increases. In reference 1, Peters and Ormiston develop matrix methods that will 
perform the Fourier operations necessary to solve linear equations, eliminating the 
need for hand calculations involving Fourier series. Several investigations 
(refs. 1-3) have used this methodology to successfully solve linear response problems. 

As opposed to linear problems in which products of two Fourier series are 
encountered, nonlinear problems retaining only second-order nonlinearities contain 
products of three series. Reference A. which considers such a problem, handles prod- 
ucts of three series by multiplying the first two series in one step and the third 
in a separate step. This presupposes that the coefficients of both of the first two 
series are known, as they are in that particular application. However, if the coef- 
ficients of only one series arc known, this method of forming the product of multiple 
Fourier series cannot be used. 

The alternative to performing successive muitipllcat ions (of two series to form 
the product of multiple series) is to derive an array analogous to the Fourier prod- 
uct matrix of reference 1. It is the objective of this Memorandum to demonstrate 
the procedure for forming such an array for any arbitrary number of series. In addi- 
tion, the application of Fourier series to performing the multiblade transformation 
(refs. 5 and 6) is described. This application uses summation relations, Fourier 
derivatives, and products of two and three series. Some of these operations can be 
found in reference 1, but a complete set has been included in this report. 


LINEAR OPERATIONS 


In reference Peters and Ormlston present the identities and linear operations 
that are reqviired to develop Fourier products, and to apply the generalised harmonic 
balance method to the solution of differential equations. For reasons of complete- 
ness and to establish a consistent notation within this report, these identities and 
linear operations are repeated below. 

Basic identities - Consider a function f(i|») with period 2n. The Fourier coef- 
ficients of are defined as follows: 
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In periodic systems having b equally spaced, identical, time-lagging components 
(e.g. , helicopter rotors, turbine compressor rotors, and wind turbines), the follow- 
ing summation relations are frequently useful. 
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in reference 1, but can be found in reference 6. They ai.^se only for periodic sys- 
tems in which b is even. Equations (5) and (6), on the otlior hand, are applicable 
to systems in which b is either even or odd. 


Phase change - To express the phase-shifted function f(<|* 4 6) as a function of 
f(l|(), let g(4») - f(^ -f 5) and express f and g as Fourier series, where N Is a 
nonnegative Integer. 

N 

fC'f') " ]C t*n + *>n •ln(n<»)) (9) 

n-o 

N 

g(*) ■ ]C tv cos(n*) + Bn sln(ni|i)] (10) 

n-o 

Writing g in terms of f and expanding, 

N 

g('f') " cosin'^. 5)] + bn airi(n(ip <■ 6)]} (11) 

n-o 


g(4') ■ ^ {(an cos(n6) + b„ sin(n5) ]cos(n^) 


n-o 

+ (-On sln(n6) + bn cos(n6) ]sin(n^) } 

The An and Bn are then 

An - Sn cos(n6) + bn sin(n6) | 

Bn " -Sn sin(n6) + bn cos(n6)/ 

Th.ls relationshl;' can be written in matrix form as 
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In the next section, a matrix operation which forms the derivative of a Fourier 
series will be discussed. That derivative matrix is just a special case of equa- 
tion (15), where the phase angle 6 equals 90”. 

Derivative- In order to apply the harmonic balance method to a system of dif- 
ferential equations, it is necessary to be able to take the derivative of a Fourier 
series. Given a serl f(<|i) and its derivative, 

N 

f (’!')“ 51 f*n cos(n4i) + sin(ni(()l (16) 

n“0 

af N 

— - ^ [-nan sin (nip) + nb„ sin (nip)] (17^ 

^ n«o 

set g(ip) equal to the derivative of £(ip). 

N 

g(ip) “ ^ [An cos (nip) + Bn sin(nip)] (18) 

n"0 


Then, 


An “ nbn 
Bn * -nan 


(19) 


Equations (19' can be written in matrix form as 


where 



( 20 ) 


( 21 ) 


By extension of the concept used to develop equation (21), the derivative matrix 
used to calculate the kth derivative of f(ip) can be obtained. 
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( 22 ) 



where 

[D*^] - [D]*^ 


Peters and Ormlston, In reference 1, discuss the case where f(^) Is nodulated 
by the harmonic function This concept can be generalized for the case where 

the Fourier coefficients of f(ij;) are also functions of ti>. 

N 

f(tj<) - 53 [an(^»)cos(ni|)) + bn(^)sin(nip) ] (23) 

n«o 


Taking the derivative of f and calling it g, 

'3a^ 3b. 
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In matrix form, where (I] is the Identity matrix. 



(26) 


(27) 


One of the uses of equation (27) will be seen during the discussion of the method of 
multiblade coordinates. 


Products of Fourier Series 

In this section, a method for forming products of any n\mbe.' of Fourier series 
using matrix algebra is described. The Fourier product is one of the most essential 
operations required to apply the generalized harmonic balance method to the solution 
of differential equations. Reference 1 discusses the technique used to multiply two 
series together, but does not attempt Co extend it Co multiple products. In order to 
form a basis of understanding of the methodology and to establish consistent notation. 
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products of two Fourier series are reviewed. Then, It Is shown how to extend the 
method to I’roducts of several series, using the product of three series as an 
example. 


Products of two Fourier series- Consider two Fourier series expansions, f(4») 
and X(ip), where the coefficients of f(i|/) are known and I and J are nonnegative 
Integers. 

I 

fC'J') “ 53 (ai cos(li|() + bi sln(l,|()] (28) 

1-0 

J 

x((|») “ 51 l*Jc cos(J't') + 5*js sin(ji|»)] (29) 

The product of f and X Is defined to be Z(^), where 

N 

Z('l') - f^nc cos(ni|») + sin(n'l')] (30) 

n»o 


where N Is a nonnegative Integer. Expanding Z In terms of f and X, 

Z('l') “ 7 ]C 23 {aiXjclcos(l + j)(j< + cos(l - j)i|»] + aiXjg[sln(l + j)i(i 
1-0 J-O 

+ sln(l - j)ip] + b£Xjg[-cos(l + j)ij) + cos(l - j)i(i] 

+ biXjc[sln(l + j)\ii + sln(l - j))**]} (31) 


The object now Is to develop a way to express equations (30) and (31) In the form 


^nc^ 

1 - [P2(f)l 

^Xjc 

2ns) 

l*js) 


(32) 


where [P 2 (f)l Is the Fourier product matrix for the product of two series. 


In order to form equation (32) the equivalent harmonics In equations (30) and 
(31) must be Identified and separated such that 1 + j - n and 1 - J - n. To this 
end, the following matrices are defined (a representing either a or b): 
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• 

-n T j <0, -n+j > I 




Note that 

the nth row of the 

Q-matrices represents the 

nth 

harmonic of Z, 

and 

the j th 

column will be multiplied by the Jth harmonic component of X. The 


harmonic 

component of a 

that 

goes 

in the nth row and Jth 

column of the 


Q-matrices is calculated according to equations (33)-(36). 

If 

the subscript of 

a 


l8 called 1, it can be seen the Q-matrlces provide a means of identifying like 
harmonics. It is important to understand how this works for the oroduct of two 
series since the same principle holds for multiple products, but is much more diffi- 
cult to produce without some sort of pattern to follow. 


The Q-matrices can now be combined to form the cosines and sines of like har- 
monics. The following substitutions are used: 


i“0 


N 


oi cos(i + j)i*' - 2^ + Q^j"(a)] ; p - i + j 


( 37 ) 


m“o 


23 Oi cos(i - j)4> “ y lnm^Qm1^“^ + Qn^(a)l ; n - 1 - j (38) 

1*0 m “0 


I 

sin(i + j)if; 
i"0 


QnJ <“) - Qnj~<«) 


n - 1 + j (39) 


I ++ _ + 

23 a± sin(i - j)ii< - " Qnj^“) 

i-o 


n - i - .1 


(40) 


where a may represent either a or b, and where [1] is a correction matrix 

■1 
2 


m - 


( 41 ) 


In equations (37)- (40), the significance of ':he pattern of signs must be stressed. 

By changing from (i + J) to (i - j) , the second sign in the C-matrix superscript 
changes from minus to plus. Because of the way in which the Q-matrices were origi- 
nally defined, the second superscript reflects how the jth harmonic of X con- 
tributes to the definition of the subscript of a. Thus, it is logical that when 
J changes sign, the second subscript should also change. In addition, the signs of 
the second matrices in equations (37) and (38), which define .ne cosines, are posi- 
tive; while in equations (39) and (40), which define the sines, the same matrices 
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art ntgatlvc. These sign paccems will be valuable when the aniltiple product 
■atrices are formed. 

All of the elements required to form the Fourier product matrix for products of 
two series have now been developed, and [^2(^)1 itaelf can now be formed. First, 
partition [P 2 (f)] auch that 



Referring to equation (31) and making the substitutions in equations (37)-(40), the 


partitions can be defined. 

[AC] - i li][tQ‘*'"(a)] + [Q‘"(a)] + [q‘*’%)] + [Q"%)1 (43) 

IA«] - \ [[Q‘^‘(a)] - tQ""(a)l - [Q'^%)] + IQ"%)]] (44) 

[B®] - i li][-[Q‘^'(b)] - lQ"‘(b)l + [Q'^^b)] + lQ"‘*(b)]] (45) 

[B*"] - -|[lQ'^‘(b)] - [Q“(b)l + lQ'*"*(b)] - [Q'%)1] (46) 


Note that since the first rows of [B^] and [A®], and the first columns of [B*] and 
[A®] are all zeros, [P 2 (f)l is singular. This difficulty can be circumvented in at 
least two different ways: an Sg may be inserted in the upper-left comer of [A*], 

or the offending rows and columns may be removed. 

Once developed, this approach to multiplying to Fourier series Is easy to use 
and highly adaptable to computer applications. When dealing with relatively short 
series, a meth^ such as this one is probably unnecessary. However, when a product 
is formed from two series each of which has many terms, or from more than two series, 
this methodology is extremely useful. 

Products of multiple Fourier series- In the preceding discussion, the steps 
required to form the Fourier product matrix were described. These steps are: 

(1) expand the product in summation notation, equation (31); (2) define the 
Q-matrices, equations (33)-(35); (3) coaO>ine the Q-matrices to define the cosine 
and sine substitutions, equations (37)-(40); and (4) partition the product matrix 
and define the partitions using steps 1 and 3. In forming the products of multiple 
Fourier series, the same steps can be followed to form arrays that can be used to 
multiply the series together using simple tensor-like operations. 

The first step of expanding the product using sumsution notation is the only one 
that requires any extensive hand calculations. This expansion is necessary in order 
to determine the signs of the cosines and sines for the substitutions chat are made 
in step 4. As an example, consider the three Fourier series, f(4»), X(v)» Y(\|.'), and 
their product Z(^), where the coefflcientr of f are known and I, J, K, and N 
are nonnegative Integers. 


fCii*) ■ S t«i C 08 (i\|i) 8 in(i<(»)] 


i-o 


X('l') - l*jc •♦■ *js 

j"0 


Y('j') - [y|tc C08(k\|<) + /i^g 8 in(kip)] 
k-a 


N 

Z('i') ■ f*nc + *08 

n«o 


( 47 ) 

( 43 ) 


( 49 ) 


( 50 ) 


Expanding Z in tenia of f, X, and Y, 

1 J K 

£ ]C H + j + '^>'1' + cos(i + j - k)t + coed - J + k)\|»] 

i“o j“o k*o 

+ cos(l - j - k)4>] + aj^xjj,yj^g(sin(i + j + k)i|< - sln(i + j - k)* 

+ 8in(i - j + k)ij; - 8in(i - j - k)4»] + a^Xjgyjj{.[sin(i + : ■*■ k)\|( 

+ sln(i + j - k)4< - 8in(i - j + k)i|/ - 3in(i - J - k)ij/] 

+ aiXjgyit5[~cos(i + j k)ij; + cos(i + J - k)i|i + cos(i - ’ + k)ij/ 

- cos(i - j - k)tjj] + b£Xj£.yitc[sin(i + j + k)4/ 8in(i + j - k)\li 
+ 8in(i - j + k)ij; + 8in(i - j - k)ii/] + biXjcyk8l"C08(i + j + k)i)' 

+ cos(l + j - k)ij< - co8(i - j + k)ii/ + cos(l - j - k)i{i] 

+ bjXjgyjjg[-cos(i + j + k)\|) - cos(i + j - k)i(i + co8(l - i + k)\|/ 

+ cos(l - J - k)ij/] + bj^Xjgyjj^g[-8in(i j + k)t(i + sln(l + < - k)i(i 

+ Bind - j + k)(J/ - Bind - j - k)i|»]} ( 51 ) 

In step 2 , where the Q-arrays are defined, the process of Identifying like hamon- 
icb Is begun. To fora the Fourier product of M series, 2 ^ Q-arrays are required. 

The first step in defining these arrays is to write all of the possible conbinationa 

of indices i, j, k, ... and set them equal to n. Then, solve each expression for 
i, which 18 the subscript of each element of the Q-arrays. For the example of the 
product of three series, the possible index combinations are: n ■ i + j k, 

n - -(1 + J k), n - i + j - k, n - -d + j - k), n - i - j + k, n - -d - J + k), 

n - i - j - k, and n - -(i - ) - k) . The eight Q-arrays (all three-dimensional) 

are then 
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i + J + k - n 


(52) 


QnjV<“> 


®n-j-k 

0 



Ofn-j-Kil 
n-j-k<0, n-J 

0<-n-j-k<I 
-n - J - k < 0 


k > I 


-(1 + j + k> 


n 


(53) 




“n-J+k 

0 


05n-j+k<l 
n-j+k<0, n-J+k>l 


1 + J - k - n 


(54) 


_ _ + “-n-j+k 

Qnjk (“) - 

0 


0<-n-j+k<I 
-n - j + k < 0, -n - J + k > ll 


-(i + j - *<) - n (55) 


(«) ■ I 


01n + j- k<I 
n+j-k<0, n+j-k>I 


i - j + k - r 


(56) 


_ “-n+j -k 

Qnjk <«) - 

0 


Oi-n + j - k<I 
-n + j - k 0, -n + J - k > l| 


-(1 - j + k) - n (57) 


+ . “n+j+k 

Qnjk (a) ‘ 

0 


05n + j+ k<l 
n + j ^ k > 1 


1 - j - k - n 


(58) 


Qnjk^<'=‘) 


“-n+j+k 

0 


Ofn + j+ kfl 
-n + j + k<0, -n + j 


}-(t - j - k) 
+ k > ll 


n 


(59) 


Note that theae eight arrays were defined by Inspection. There was no need to %n:lte 
out each one to check Its contents, because they all follow a specific pattern which 
defines their use. 


Step three completes the Identification of like harmonics by combining the 
Q-matrlcea to obtain the cosines and sines of each harmonic. Again, this can be 
done by Inspection. Starting with cos((i 4 - j + k + ...)i(/j, the Q-arrays that rep- 
resent n- (i+j +k+ ...) are added together and premultiplied by a correction 
matrix which multiplies all of the n ■ 0 terms by one-half. Similarly, the other 
cosine functions are obtained by adding the appropriate Q-arrays and premultip lying 
by the correction matrix. The sine functions are defined by stibtracting the 
Q-arrays from one another and no correction is i<;>;uired. Returning to the exaa^rle 
of the product of three series, the cosine and sine expressions are 

1 N _ 

5^ cosd + j + k)ij- - £ InmlQnjk + Qmjk J n - i + j 4- k (60) 
1-0 m-e 
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i M 

^ co8(i + J - k)i|; - ^ 
l«o m*o 

I N 

a, cos(i - j + k)4» •» ^ 

l-*o o«o 

I N 

2^ cos(i - j - k)\(i - ^ 
l-»o m»o 


InalQnJk ■*■ ^mjk * 


n - i + J - k (61) 
n - i - j + k (62) 
n • i - j - k (63) 


2^ aj cin(i -v*- 
1 *0 

1 

23 aj Bind + 
i“0 

I 

^ 8in(i - 
i-o 


+ k).*. - QnjV^“> 

- - 3njV<“> 

+ k)4- - Qnjk~^“) 


Qnjk 


QnjV<“' 


QnjV<®> 


n - i + J + k (64) 
n - i + j - k (65) 
n - i - j + k (66) 


23 “i Bind - j - k)ii) « Qnik^^®^ " ®nik^^“’ » n » i - j k (67) 

l-o •' 

whero o may represent either a or b, and where I is the correction matrix 
de '^led in equation (41) 


The final step in forming the product of multiple Fourier series ia to parti- 
tion the product array [Pm(f)] and define tne partitions. Within tPi.j(f)] the cosine 
and sine coefficients are grouped together in blocks which are arranged in a checker- 
board fashion throughout the array. If the partitions are chosen such that each one 
Contains only cosine or sine coefficients of f(i|;), there will be 2^ partitions. 
Writing the relationship among Z(\^>)> ff’l'lt X(4») • and Y(i^) for the example of a 
product of three series in pseudo-matrix notation. 


“nc 


no' 


iPj(f)] 


^jc 
’'jsJJ 


[Ykc 

Ivks 


Partitioning [P,(f)1. 


I J“0 j"0 ^ 


z_„ ■ 5^ I 52 *4 )yt. 

k-e I J-o J8 -^kc 


J 

jTo 


(^^8) 


(69) 


(70) 


Using equation (51) to define the partitions and substituting equations (60)-(67), 
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*i>Jk ■ 4 ^ ‘Vjk ‘*«lk*‘*' '*' **ijk**‘* ■*■ “Cik **' 


* \3k’<‘> ^ «,Jk <•> * <Wk 0)1 


,C8 _ 1 r^+ - -k 


4 - Qnik <«> - Qn^v <“> + Q”7 i, (a) + Q„, "(a) - Q‘ “(a) 


njk 4 '^njk 


- <jVo) 4 q;;/u)i 


\jk * 4 ^ ^nmf-Q^jk <®> \jk ^jk'^®^ 








B„j\ - i r 




®njk " 4 ? ^nm[-0;jk ■*■ 




“njk - i f-Ck"' 
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Having obtained the partitions of the product array, they can be substituted into 
equations (69) and (70), and (P8(f)] can be formed. Again, care must be taken in 
using [Ps(f)] because each "slice” of that array (j * 1, 2, ... J or 
k - 1, 2, ... K) is a singular matrix, as described in the section on the Fourier 
product matrix. As in that case, the problem can be eliminated quite easily. 

Forming the product array for products of M Fourier series can be quite 
simple. However, it is obvious that the more series that are multiplied together, 
the more cumbersome this process becomes. Fortunately, once the product array has 
been derived for M series, it need never be done again. 


Multiblade Coordinate Transformation 


The purpose of the multiblade transformation is to transform equations from a 
coordinate system which rotates at a constant angular velocity to a nonrotating 
coordinate system. In reference 5 the mathematical definition of multiblade coordi- 
nates is presented, but not thoroughly explained. Reference 6 contains a more com- 
plete development of the mathematics of the transformation and its use. Since the 
multiblade transformation is really only a modified Fourier series, the operations 
that were defined in the preceding sections should be applicable after making some 
minor modifications. It will be seen that using certain Fourier operations makes 
the transformation simple to perform for any number of blades, with no ad hoc addi- 
tions or corrections. 


Consider the linear (or linearized), second-order, differential equations for a 
dynamic system having b identical, equally spaced, time-lagging components. The 
kth component has a degree of freedom qj^ which is referenced to the rotating 
coordinate system. There is also a degree of freedom p which is referenced to the 
nonrotating system and is common to all b components. In the rotating system, the 
equations of motion for qj^ and p are respectively 


”qqWk>’k * ‘^qq<*k>’k %q«k>’k * "qp<V'’ ‘^qp‘*k>'’ 

”pq‘*k>’k Sq‘*k>’k “pq^^k’^K * "pp'^’P * *=pp<*k>P 

where the coefficients of the equations are periodic in 


+ KqpWk'P - Fq<1'k> 

(79) 

+ ■ Fp(»k> 

(80) 


4'k “ 2ir(k - l)/b. 


If these equations were to be put into the form of a variational statement, 
equation (79) would be multiplied by 6qj^ and equation (80) would be multiplied by 
6p. The resulting variational statement Is 


t 


{[M q, + C q, + K q, + M p + C p f K p - F ]6q. 
qq^k qq^k qq^k qp*^ qp*^ qp*^ q ^k 


+ M q. + C q, + K o, + M p+C p+K p-F ]6p}dt - 0 (81) 

pq k pq^k pq k pp^ pp*^ pp*^ p 


To transform equation (81) from the rotating system into the nonrotating system, only 
6qk, qj^, and Its derivatives need to be transformed. The degree-of- freedom p is 
already referenced to the nonrotating coordinate system. Since derivatives of the 
coordinate qj^ would be difficult to manipulate, the dummy variables and 
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trill b« substituted for end qj^. respectively. Applying the definition of the 
multiblsde coordinates (ref. 5), 

N 

- Oo + Od(-l)^ + 53 [oinc cos(n*it) (»„, sin(n<|>i^)] (82) 

n-i 

N 

Bjt - Bo + E tBnc cos(n^jj) + B„s sin(n*^)] (83) 

n»i 

N 

<lk " «lo + ]C ^*lnc co»(n'l'k) ‘Ins »in(nH»k)l 

n-i 

N 

fiqk “ *^*10 D f^‘»nc <io®<n'l'k^ + fiqns Pin(n^)] (85) 

n*i 

Since the coefficients of equation (81) are periodic functions (i.e. , Fourier 
series), substituting equations (82)-(85) into that equation results in terms that 
are products of modified Fourier series. The terms ^qq^k^qk* ^qq®k^qk* 

Kqqqk^qk products of three series, Kppp6p, Cppp6p, Kpppdp, and FpSp are 

single series, and all of the others are products of two series. 

To form the products required in equation (81), the same basic methodology 
described in the preceding section can be applied. The difference lies in the 
treatment of the (-l)I^ terms. Since the locations in the arrays allotted to the 
sin( 04 >) terms are always filled with seros, those locations may be used for the 
(-1)^ terms. This also preserves the synnetry of the arrays. 


The formation of the modified Fourier products in equation (81) performs the 
transformation from the votating to the nonrotating coordinate system. To complete 
the process, the dummy v ^iables a]^ and 8)^ need to be eliminated. Using the defi- 
nition of the derivative of a series having periodic coefficients 


\ 


0 

®1C 

Ojc 


qo 

qic 

q?c 


“nc I 
“d 
'is 
‘28 








qnc I 
‘Id 
qis 
qis 


( 86 ) 


*ns , 


1*108 / 


'Po ' 


% 

Pic 



P2C 

• 

• 


^2 

• 

• 

• 

»nc 

■ • d » ml 

• 

*lnc 

0d 


*ld 

eis 


qis 

Pzs 

• 

• 


^28 

• 

• 

• 

‘ ®ns‘ 


• 

‘‘^ns * 


Once equations (86) and (87) are substituted Into the transformed equations, the 
transformation is essentially complete, except for sorting the periodic terms. 

Thus far, the number of rotating components, b, has not been needed for the 
transformation. However, it is of great importance in determining the number of 
multiblade coordinates that are needed, as well as determining the harmonic contri- 
butions in the nonrotating system. The number of multiblade coordinates that are 
required is equal to the number of components, where N * (b - 2)/2 for b even 
and N ■ (b - l)/2 for b odd. If b is odd, the differential collective mode 
does not exist (ref. 5). To determine the harmonic contributions in the nonrotating 
coordinate system, equation (81) is summed over the b components and the identi- 
ties in equations (5)-(8) are applied. The transformation is now complete, and the 
variational statement may be broken down into differential equations. 

In applying this method for using the multiblade transform, there are shortcuts 
that can be taken. There is no practical reason to convert the differential equa- 
tions to a variational statement and then convert back. This was done to show the 
logic behind the ope''ations that were performed. Similarly, there is no need for the 
dummy variables and which were included only for convenience. The sorting 
of harmonics can be done by inspection. If b is odd, only those harmonics that are 
integer multiples of b are retained. If b is even, harmonics that are integer 
multiples of b are retained for all terms not involving (-1)^. In those terms 
where (-l)l^ appears, only those harmonics which are odd integer multiples of b/2 
are retained. 


CONCLUDING REMARKS 


In this Memorandum, the development of Fourier operations which was begun in 
reference 1 is extended. Summation relations for expressions involving (-1)^ are 
defined, the definition of the derivatives of Fourier series are generalized to 
include series having coefficients that are functions of ij>, and a method for per- 
forming the symbolic multiplication of any number of Fourier series is developed. 

In addition, the operations are applied to the multiblade coordinate transformation. 
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While the harmonic balance method may not be suitable for solving all types 
of response problems, its elegance and understandabllity are desirable features. 

Its major disadvantages are the large number of algebraic equations that must be 
solved and the necessity for forming products and derivatives of Fourier series. 
Because of the work performed in reference 1 and the extensions included in this 
report, the second disadvantage becomes unimportant since those tasks can be per- 
formed by the computer. 

Other applications of Fourier operations are possible. In this report, the 
example of using these operations to perform a coordinate transformation was devel- 
oped. The result of this development was a procedure which has general applicabil- 
ity, is understandable, and is suitable for automation. Many o*-her applications for 
these operations may be found. 
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